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Abstract 

We report an observation of charge conservation breaking in a model study of electronic current noise of transport through a 
dissipative double quantum dot within generalized master equation formalism. We study the current noise through a double 
quantum dot coupled to two electronic leads in the high bias limit and a dissipative heat bath in the weak coupling limit. Our 
calculations are based on the solution of a Markovian generalized master equation. Zero- frequency component of the current 
noise calculated within the system, i.e., between the two dots, via the quantum regression theorem exhibits unphysical 
negative values. On the other hand, current noise calculated for currents between the dots and the leads by the counting 
variable approach shows no anomalies and seems physically plausible. We inquire into the origin of this discrepancy between 
two nominally equivalent approaches and show that it stems from the simultaneous presence of the two types of baths, 
i.e., the electronic leads and the dissipative bosonic bath. This finding raises interesting questions concerning conceptual 
foundations of the theory describing multiple-baths open quantum systems widely encountered in nanoscience. 
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1. Introduction 

Recent advances in technology, fabrication, and 
measurement of mesoscopic semiconductor devices 
with ever-decreasing dimensions of achievable nanos- 
tructures stimulate also theoretical studies of physical 
phenomena determining their properties. Questions 
of prime interest concern their possible quantum be- 
havior and the quantum-classical crossover due to 
interaction with surrounding environment causing dis- 
sipation, relaxation, and dephasing [1-3]. One of the 
conceptually simplest and experimentally achieved 
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systems is the double quantum dot (DQD) which can 
be tuned into a regime where it is effectively described 
as a tunable two-level system for the electronic en- 
ergy states [4]. These energy states can be tuned, for 
instance, by means of an external gate voltage. The 
interest in such devices stems also from the attractive 
possibility to utilize them as potential q-bits. 

In the setup consisting of a DQD, the role of the elec- 
tronic coherence between the two spatially separated 
electronic states corresponding to the respective dots 
is of central importance. The DQD device loses coher- 
ence due to the coupling to noisy environment (e.g. 
noise in the gate voltages and unavoidable interaction 
with phonons in the substrate). Moreover, energy can 
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be exchanged with bosonic degrees of freedom which 
gives rise to transitions between states of nonequal en- 
ergy and, thus, relaxation. The dissipative dynamics of 
two level systems (spin-boson problem) have been sub- 
ject of study for many years [1,5]. DQD setup brings 
along a new twist to this standard problem in the fact 
that the DQD is electrically contacted by leads and 
charge transport through the DQD occurs [4]. This is 
a new feature of the old problem adding complexity 
to the methods of its solution. On the other hand, for 
most practically interesting setups the dissipative cou- 
pling may be considered as rather weak which should 
bring about important simplifications for its solution. 
Yet, as we will demonstrate in this work, the simulta- 
neous presence of the two different kinds of baths (elec- 
tronic leads and standard dissipative bosonic bath a la 
Caldcira-Leggett [6]) causes serious conceptual prob- 
lems within the simplest possible formalism of Marko- 
vian generalized master equations (GME) when it is 
applied to the weakly dissipative DQD problem. 

It has been suggested by Aguado and Brandes [7] 
that considering the electronic current noise in the 
DQD devices may be a useful tool for characterizing 
their dissipative properties going beyond the infor- 
mation available from the stationary current charac- 
teristics only. In this work we adopt their model and 
perform an exhaustive comparative study of the eval- 
uation of the current noise based on GME approach 
used in previous studies [7,8]. We surprisingly find an 
internal inconsistency of the formalism which breaches 
fundamental physical law of charge conservation. Two 
nominally equivalent approaches for the calculation 
of the zero-frequency component of the current noise 
spectrum, namely the quantum regression theorem 
(QRT) and counting variable approach with the Mac- 
Donald formula, show mutual discrepancy. This fact 
puts under question the status of the results obtained 
by these methods. 

The structure of the paper is the following. In Sec. 2 
we introduce our model and in Sec. 3 we describe the 
method of its solution via an approximate Markovian 
GME. The following Sec. 4 describes the evaluation of 
the current noise within the used GME formalism while 
Sec. 5 discusses the general aspects of the charge con- 
servation and its breaking within the GME formalism. 
Sec. 6 presents an overview of results for our studied 
model of a dissipative DQD. We discuss the obtained 
results and their implications together with resulting 
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Fig. 1. Schematic depiction of the studied DQD system. 

open problems and outlook in the concluding Sec. 7. 



2. Model 

The double quantum dot device [4] in Fig. 1 is de- 
scribed as two electronic levels (corresponding to par- 
ticular single-electron levels within the transport win- 
dow of the left and right dot, respectively) de-aligned 
by an energy difference e with a coherent interdot tun- 
nel coupling f2. The system is assumed to be in the 
regime of strong Coulomb blockade so that only three 
states play a role: no extra electron |0) on the whole 
DQD system, one extra electron on the left dot |L) and 
one extra electron on the right dot |R), i.e., we exclude 
multiple occupancies of the DQD. This can be achieved 
by a suitable gating, when a very high charging en- 
ergy prohibits am addition of more than one electron. 
We also consider spinless electrons. Hamiltonian of the 
DQD device then reads 

Hs = i£(|L)(L| - |R)(R|) + n(|L)(R| + |R)(L|). (1) 

The so called device bias e can be tuned by gat- 
ing. The term proportional to Q. enables the tun- 
nel current through the device. The eigenvalues of 

the isolated system Hamiltonian are _Ei,2 = 

with A = V AQ? -\- and the corresponding cigen- 

vectors read |1) = ^/^|L) + y^lR), |2) = 

-\/i?|L> + yii^iR)- 

The double quantum dot is coupled to two leads with 
a high bias applied between them. The bias is smaller 
than the charging energy but otherwise it is the largest 
cucrgy scale in the model. We assume that the non- 
interacting leads are coupled via standard tunneling 
terms 
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He + Hcs = ^ Eki^clj^Cki, + ^ Eknclj^Ckn 

k k 

+ E^*44l|0)(L1 + |L)(0|c,l) (2) 

k 

+ T.Vkn{cU0){R\ + \R}{0\ckH). 

k 

The leads are held at respective electrochemical poten- 
tials /iL and fin whose difference gives the bias. We as- 
sume that yUL — > 00 and /xr —^ — oo. The tunneling den- 
sities of states Ta{e) = ^J2k |Vfca|^(5(e - Eka) ,a = 
L, R are assumed energy independent (wide-band limit 
or first Markov approximation, compare Ref. [2]) and 
equal Fl = Fr = F. Both the high bias and wide-band 
limits are necessary for the applicability of the Markov 
approximation later on. 

Finally, we introduce a generic dissipative heat bath 
a la Caldeira-Leggett consisting of an infinite set of 
harmonic oscillators [1,6] which are linearly coupled to 
the left-right population difference of the DQD [7] 

Hb + Hbs = ^iio-lo-j + l) 

' (3) 
+ Y.Ci{a]+a,){\h){h\-\K)m). 

i 

The heat bath is fully characterized by its spectral den- 
sity 

J{u>)=2Y,\C,f5{Lo^uj,), (4) 

3 

which wo take in the Ohmic form J{ijj) = 2h?^ijj/-K ■ 
exp (— w/ojc). The parameter 7 gives the strength of the 
dissipation and Wc is a high energy cut-off frequency 
[1,6]. 



3. Generalized master equation 

3.1. Liouville space 

In order to conveniently manipulate with density 
operators, we define Liouville space [8]. The Liouville 
space is a linear space spanned over operators acting 
on the original Hilbert space assigned to the system. 
Its basis |n, n')) is constructed from a basis jn) of the 
Hilbert space as \n,n')) = \n){n'\. A general operator 
in the Hilbert space 

A= '^A„„'\n){n'\ 

7i,n' 



corresponds to the vector 

\A}) = ^A„„,\n,n'}} 

in the Liouville space. For example, a density operator 
p can be written as p = Yln,n' Pnn'\n, n')) with p„„/ = 
(njpjn'). We define the scalar product on the Liouville 
space by {{A\B)) = Tr{A''B}. In order to avoid confu- 
sion, linear operators acting in the Liouville space are 
called supcroporators and in the following will be de- 
noted by calligraphic symbols. The vectors of the Liou- 
ville space in the brarket notation will be distinguished 
by double brackets. The matrix representation of su- 
peroperators then follows from the previous 

^ = ^ A„n',mm'\n,n')){{m,m'\. 

n,rL' 

3.2. Liouvillean 

For the description of a dissipative system we dis- 
tinguish between the system itself (electronic states 
of the dots) and the reservoirs (heat bath and leads). 
Our task now is to get a closed evolution equation for 
the reduced density operator which is the system part 
only of the total density operator. To this end we per- 
form the standard projection onto the system assum- 
ing weak coupling to the reservoirs and consequently 
using the Markov approximation [2] . Due to the weak 
coupling the effects of the two baths are additive. 

It should bo noted that within our assumptions on 
the leads (wide band limit and high bias) an equiva- 
lent result for the effect of leads can be obtained with- 
out the weak coupling assumption, i.e., for arbitrary F, 
as shown by Gurvitz and Prager [9]. It turns out that 
these assumptions correspond exactly to the so called 
singular coupling limit in the mathematical literature 
(see, e.g., Ref. [10]) which also leads to a Markovian 
dissipative evolution of the system. The additivity of 
the two bath is then, however, only a heuristic assump- 
tion which may break down for large enough F. The 
range of parameters we consider is restricted to rather 
small F so that the potential differences between the 
two possible approaches should be safely negligible and 
the result of the projection leads to the following GME 

^ = Cp{t) = Csp{t) + Cnpit) + Ccp{t) , (5) 
with 
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Csp{t) = 



--[/?s,p(i)], 



At 



X TrB{[HBS,[HB3i-T), p{t) pb]]}, 



X Tic{[Hcs, [Hcs{-t), p{t) ® pc]]}. 



The part £s describes the free evolution of the sys- 
tem while jCb, jCc determine the dissipative influence 
of the generic heat bath and the electronic leads, re- 
spectively. It turns out that the off-diagonal elements 
Pok, Pko with = 1, 2 of the reduced density matrix are 
decoupled from the rest of the system, i.e., their evo- 
lution does not enter expressions for the other matrix 
elements and vice versa (compare with Ref. [8]). There- 
fore, the subspace {|0a)), |qO))} {a — L,R) can be pro- 
jected out leaving us with the relevant Liouville sub- 
space with the basis {|00)), |LL)), |RR)), |RL)), |LR))}. 
In this basis the above parts of the total Liouvillean 
are described by the following matrices [7,11] 



4. Current noise 

In this section we briefly introduce the current noise 
[12,13] and methods of its evaluation within the GME 
framework. More detailed account of these issues can 
be found in Sec. Ill of Ref. [8]. 

Equations of motion for the operators of the occu- 
pation of the loft dot riL = |L)(L| and the right dot 
nR = |R)(R| read 

d ic 
e-r-nj, = ~-r[ni,,H] = - Inh, (H) 
at h 

d ic 
e ^riR = - — [nn, H] = 7rl - /or- (12) 

On the right side of the equations, we identify charge 
current operators across the different junctions: /lo = 
— ^[riL, -Has] is the operator of the current between 
the left load and the left dot, 7rl — ■^[nL,-ffs] = 
— -^[fiR, Hs] is the operator of the current between the 
dots, and /or = j; [rm, Hcs] is the operator of the cur- 
rent between the right dot and the right lead. Explic- 
itly, they read 



= f E^-^i^I^lIOXLI - |L)(0|c.l), (13) 



Cs = T 



Cb = 



Cc = 
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(6) 



(7) 



(8) 



TT n 

7± = -T^ 



^|j(A/R)[Jcoth(i/3A)±l], (9) 

(10) 



/RL = |n(|L)(R|-|R)(L|), 



^01^ = I E (I <0| '^^R - 4r 1 0) {R| ) . 



(14) 
(15) 



Since the commutators with the bath operators are 
zero [nL, Hbs] = ['^R, Hbs\ = 0, the heat bath gives no 
explicit contribution to the current operators. 

The current operator /rl is obviously a system op- 
erator, i.e., it acts as unity on the degrees of freedom 
of the leads and the heat bath. However, this is not the 
case for the two operators of current between the dots 
and the leads /lo and /or. 

Next, we define the current autocorrelation function 



CA(r)= lim [l{{U{t + r),lA{t)}) 



(16) 



-(/A(t + r))(/^(t)>], 

with A = LO, RL, OR. Due to the stationary limit (i — > 
oo) the autocorrelation function is symmetric Ca{t) = 
Ca{—t). We define the current noise spectrum as 



Sa{oj) ■ 



J — C 



drCA(T)e'"". 



(17) 



The current noise spectrum is non-negative as can be 
shown by using the Lehmann representation. 



4 



Now we need to express the current noise spec- 
trum in terms of the quantities involved in the GME 
(5). We denote the stationary reduced density matrix 
limt^oop(t) = Pstat = |0)). It satisfies Cpatax = 0, 
hence it is the zero-eigenvalue (right) eigenstate of the 
Liouvillean. Since the Liouvillean is not Hermitian, left 
zero-eigenvalue eigenstate denoted by ((0| is not just 
the Hermitian conjugate of the right zero-eigenvalue 
eigenstate |0)). However, one can see that ((0| = 1, 
because for an arbitrary system operator A 

((0|£|A)) = Trs(l£yl) = Tis{CA) = 

due to normalization of the reduced density matrix. 
Now, we define the projector on the kernel of the Li- 
ouvillean V = |0))((0| and its orthogonal complement 
Q= l — V. With help of Q the well-defined superoper- 
ator TZ = QjC~^Q represents the pseudoinverse of the 
Liouvillean ( "inverse on the regular subspace Q" ) . 

With that we have all necessary ingredients for ex- 
pressing the current noise. It can be done in two dif- 
ferent ways for the two types of junctions. The RL- 
junction lies in the system and, thus, the quantum re- 
gression theorem (see Ref. [2], Sec. 5.2) can be used 
to calculate the correlation function C7rl(t). The final 
formula for the zero-firequency current noise reads [8] 



S'rl(O) = -2e^((0|2:RL7eXRL|0)), 



(18) 



where we have introduced the current superoperator 
1 



2^RLp = ^Url,p} 



(19) 



in terms of which the stationary current is given as 
(/rl) = Trs(/RLP.tat) = e((6l2:RLlO)). 

For the outer junctions (between the dots and the 
leads) the QRT cannot be used, because the current 
operators 7lo and /or involve the lead operators. How- 
ever, n-resolved form of the generalized master equa- 
tion and the MacDonald formula [13,14] enables us to 
calculate the zero-frequency noise also for these junc- 
tions, see details in Ref. [8], Sec. III. For the stationary 
mean current through the OR-junction we get (/or) = 
eTrs(2:oRPstat) = e((0|ToR|0)) with 

IoRP = r|0)(R|p|R)(0|. (20) 

The final result for the zero-frequency current noise 
reads 



S'or(O) = e2{{0]2:oR - 22:oR7^ToR|' 



(21) 



and analogously for the LO-junction with TloP 

r|L)(0|p|0){L|. 



5. Charge conservation issue 

The equations of motion (11) and (12) for the dot 
occupation operators (charge conservation conditions) 
imply that the stationary mean current and the zero 
frequency noise are independent of the measurement 
position along the circuit [8] 

(/lo) = (/rl) = (/or) , 5'lo(0) = S'rl(O) = 5or(0). 

Let us now focus on the reformulation of the charge 
conservation condition in the superoperator language 
and evaluate, e.g., the commutator [JVl,/^] with the 
superoperator of occupation of the left dot defined by 
JVlP = i{nL,p} in analogy with other superoperators 
corresponding to system operators such as, e.g., Trl. 
Its matrix representation in the relevant 5-dimensional 
Liouvillc subspace reads A/l = diag(0, 1, 0, 1/2, 1/2). 
Then with the help of Eqs. (6), (7), and (8) we arrive at 



(22) 



with the "anomalous current superoperator" Ia = 
— [A/ljJCb]- Since [nL,//Bs] = 0, we would expect that 
also [jVl, jCb] is zero, therefore the final Liouville space 
analogy to the charge conservation condition (11) 
should read [A^l, C] = Tlo — Irl. Indeed, if this condi- 
tion were satisfied one could easily explicitly show the 
equivalence of Eqs. (18) and (21) using the method of 
Ref. [8], Sec. HID. Unfortunately, as one can readily 
evaluate 

/O 0\ 


, (23) 

7+ 7- 
\0 7-H 7- 0/ 

so that the anomalous current superoperator is not 
identically zero. Analogously, for the other dot we ob- 
tain the result 



lAfR,C] ^Irl +Ia-Io 



(24) 



with the opposite contribution form the heat bath com- 
pared to the left dot, i.e., [A/'r,-Cb] = — [A/lj-Cb] =^a- 
Thus, we face the fact that the charge conservation 
condition contains nonzero anomalous terms. The pos- 
sible consequences may be that the mean current or 
the zero- frequency current noise are no longer equal for 
all pairs of junctions. First, since ((0| corresponds to 
the unity operator and, thus, is equal to (1, 1, 1, 0, 0), 
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we see that ((0|Ta|0)) = and the mean current is con- 
served along the whole circuit. Also, we notice that 
[A/l +A/'r, jCb] = implies there is no problem with the 
charge conservation between the outer junctions. How- 
ever, the zero- frequency noise eventually does show dis- 
crepancy between the outer junctions and the inner 
one, in particular 

SrlCO) - Slo(O) = e^((6|TLo7^2:A|0)) = e'rae 
^ h^r^+A' + 2n^A(7+ + 7-) - 2£^7+(f 7- - 



[ire2 -I- 3n2A + iR2rA2 - Qs{2-y+ + 7-)] 



(25) 



where A = 7p + ^F. We will analyze features of this 
discrepancy in detail in the next section. 

Before that, let us return to the point where we have 
assumed that [A/lj-Cb] was identically zero. This was 
a reasonable assumption, because the dot occupation 
operator til commutes with the heat bath-system in- 
teraction Hamiltonian Hbs and thus the bath vari- 
ables do not enter explicitly the current operators. It 
can be shown that [A/l , jC] indeed does not depend on 
the bath variables if we operate on the whole Liouville 
space (system -|- bath) before projecting onto the sys- 
tem and introducing the weak coupling and Markovian 
limits [15]. However, after the projection we arrive at 

[A/'l,£b]p = 

1 r°° 

= — y dTTVB{{[/?BS,nL],[/?Bs(-r),p®pB]}} 

1 f°° 

+ drTrB{[i?BS,{[i?Bs(-r),nL],p®PB}]}. 

(26) 

Apparently, the first term is equal to zero due to 
[nh,HBs] = as expected, while the second term 
yields the nonzero anomalous current. 



6. Results 

In this section we present various aspects of the pre- 
viously found charge conservation breaking in the dis- 
sipative DQD model. In the following we will conven- 
tionally represent the noise by a dimensionless quan- 
tity, the Fano factor F = S{0)/e{I) [12]. 



6.1. Zero dissipation 

In the case without the heat bath (7 = 0), the charge 
conservation condition is satisfied. For the mean cur- 
rent [16] and the Fano factor [14] we obtain 

n2 



(7) = eV 



£2-H3n2 + (l/,,r)2 ' 

(O^ + 2(iftr)2) 



(e2+3Q2+(l;jr)2)^ 

These results arc illustrated iu Fig. 2. The mean current 
vs. bias £ has the Lorentzian shape with the half-width 
^/3QFnWJi)^ and maximum at e = 0. The Fano 
factor has the dip at e = where quantum colioroncc 
strongly suppresses the noise. For large |e| > the 
mean current becomes very small and thus electrons 
tunnel very sparsely and consequently the tunneling 
events are uncorrelated which corresponds to a Poisson 
process with the value of the Fano factor F — » 1. 

6.2. General case 

When dissipative heat bath comes into play (7 > 0), 
the transport is strongly affected by the possibility of 
exchanging energy with the heat bath [7] as it is illus- 
trated in Fig. 2. The shape of the mean current curve is 
no longer Lorentzian but exhibits an asymmetry. With 
increasing temperature the peak becomes broader and 
more symmetric. Analytically, we obtain 

era (QA - 7+e) 



(I) 



\Te'^ + 3fi2A + iR2rA2 - ^£(27+ -|- 7. 



(27) 

In the right plot of Fig. 2 Fano factors -Flo ~ For (solid 
lines) and Frl (dotted lines) are plotted. The differ- 
ence between the Fano factors obtained by different ap- 
proaches is significant. Interestingly, at £ = the Fano 
factors have the same value Flo = -Frl as follows from 
relation (25). 

For £ > spontaneous emission occurs even at very 
low temperatures and the noise is reduced well bel- 
low the Poisson limit. Larger couplings 7 lead to very 
asymmetric Famo factor. At finite temperatures, ab- 
sorption of energy quanta from the bath is also pos- 
sible and the Famo factor for £ < is reduced bellow 
the Poisson limit too. With increasing temperature the 
effect of the emission and the absorption is growing, 
except the point £ = where both the mean current 
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Fig. 2. Mean current and Fano factor vs. device bias £ 
for different values of damping coefficient 7. Solid lines 
show the mean current and the Fano factor at the outer 
junctions LO, OR, dotted lines at the inner junction RL. 
Parameters: Q = 5meV, F = 0.1/ftmeV, /3 = 0.1 meV"!. 





[meV] 



Fig. 3. Fano factor vs. device bias e for F ^ 0. Solid 
line shows i^LO = ^ORi dotted line i^RL- Parameters: 
n = 5meV, (3 = 0.1 meV'^. 



6.3. Limit V ~* 



and the Fano factor are temperature independent. It 
appears that the MacDonald formula yields physically 
plausible results for Flo and Fob., whereas Ful given 
by the quantum regression theorem behaves patholog- 
ically with unphysical negative values and non-Poisson 
limit for e — > 00. For sufficiently strong coupling (7 ~ 
10~^) Frl drops to negative values in the e > re- 
gion and for sufficiently high temperature (T ~ 200 K) 
also in the e < region. Analyzing the expression (25) 
we find that for e ^ 00 the noise difference AS = 
Srl(O) - S'lo(O) ~ 1/e. From the relation (27) for the 
mean current in the same limit it follows (7) ~ 1/e. 
Therefore their ratio yielding the difference of Fano fac- 
tors A_F = Frl — Jlo does not go to zero for e 00 
as it should. Nevertheless, despite of the fact that T^rl 
behaves manifestly wrong, we do not have any valid 
proof yet, that the MacDonald formula gives a better 
and more reliable results for Jlo and For. 

In the following subsections we will investigate how 
both the MacDonald formula and the quantum regres- 
sion theorem approach behave in several approxima- 
tions or limit cases. It will answer whether the Mac- 
Donald formula gives physically acceptable results and 
will show more pathologies of the quantum regression 
theorem results. 



Limit r ^ could be potentially interesting - 
analogously with the dissipationless limit 7^0 the 
charge conservation could possibly be recovered. The 
mean current and the zero-frequency noise go to zero 
in this limit, however their ratio, the Fano factor, does 
not [15]. The difference between the Fano factors of 
the outer and inner junctions reads 

Flo - Frl = 

2e^ (coth^ (i/3A) - 1) 
^ £2 + 4eAcoth (i/3A) 4- 3A2 coth^ (|/3A) ' 

Results are illustrated in Fig. 3. We note quite an in- 
teresting phenomenon that the Fano factor does not 
depend on the heat bath spectral density and, thus, it 
is not inffuenced by the strength 7 of the dissipation. 
Nevertheless, all the anomalies survive. The Fano fac- 
tor T^RL can be negative for certain e and temperature 
high enough. Since the Fano factors for the three junc- 
tions are not equal, the charge conservation condition is 
not fulfilled. Both Flo and Frl become 1 for e — > —cx, 
1/2 for e ^ cxj and F = 5/9 for e = 0. The striking 
difference between the Fano factors for the inner and 
outer junctions is that Foh has no maxima or minima 
and just smoothly decreases from 1 to 1/2, whereas 
Frl has two minima, one for e < and the other for 
e > 0. 
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E |meV] 



identical results for both the mean current as well as 
the current noise. Thus, it seems that all problems are 
fixed. However, the physical content of our results has 
undergone great changes. We have given reasons why 
the mean current and the Fano factor curves should 
have emission-absorption asymmetry, but we obtain 
absolutely symmetric curves - the mean current gained 
the Lorentzian shape and the Fano factor has no sup- 
pression for e > due to the emission process. Because 
of these reasons we must reject the rotating wave ap- 
proximation on the physical grounds. 

6.5. Pauli master equation 



Fig. 4. Mean current and Fano factor vs. device bias e for 
different values of damping coefficient 7. Solid lines corre- 
spond to the outer junctions without the RWA while dotted 
lines correspond to all junctions within RWA. Parameters; 
r = 0.1/;imeV, Q = 5meV, f3 = O.lmeV-i. 

6.4. Rotating wave approximation 



The rotating wave approximation (RWA) is under- 
stood as neglecting terms of the form pa\ apa and 
keeping terms of the form pa, apa^ in the generalized 
master equation dp/ dt — Cp [2]. In our language, we 



|L)(R1 



= |R){L| and, therefore, |L)(L| = 
Application of the RWA to (7) 



o-V, |R){R| = 
leaves us with 



^±1 



The anomalous current (23) is now identically zero, 
since [Ni,,C^^] = [Nn,C^^] = 0, and therefore 
the charge conservation is restored. The zero-frequency 
noise and the Fano factor for all junctions are equal. 
The mean current is 

eFfi^A 



/O 





\ 























-7p 





^0 





-7p/ 



{!) = 



and for the Fano factor we get 

[A^ (n^ + ift"r(7p + F)) + \ij,Ve^] 



1 



These results are illustrated in Fig. 4. 

We can see that negative Fano factor does not ap- 
pear in this approach. Furthermore, all junctions give 



Pauli master equation is yet another modification 
of the original GME which does not lead to formal 
inconsistencies. When the coupling to the leads F is 
small enough, which is the case in the regime we con- 
sider here, we can neglect the off-diagonal elements 
pi2, P21 of the density matrix in the system eigenba- 
sis [10] . If we transform the Liouvillian of the original 
GME in the left-right basis (5) into the system eigen- 
basis {|0), 1 2)}, we can get the Liouvillian (or rate 
matrix) of the Pauli master equation immediately just 
by restricting ourselves to the subspace spanned by 
Poo, Pii, P22 which leads to 

/—Fin — F20 Foi 

-Foi - -Ji 
7i 



-Fio 

Fio 
F20 



F02 
7t 
-F02 — 



7T/ 



where 

Fio ~ F02 



A + e 



2A ' 

7T = ^S^(AA)' 



F20 
1 



2A 



-0A ' 



h? A2 " e«^ - 1' 

Identical results can be derived directly from a rate 
equation approach to the occupations of the eigenstates 
only with the rates determined by the Fermi golden 
rule. The current superoperators in the new Liouville 
subspace are given by 



1m) — 




2oR — 



/ Foi F02 

yO 



while no analogy of the superoperator Xrl exists on 
the chosen subspace. 
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For the mean current and the Fano factor between 
the dot and the lead, we get regardless of the choice of 
the left or right junction 

^ ' ire2 + 3n2A-fie(27++7_) 

and 

F = 1 - 20 X 

20^ A' - n'£(7p + |r)(37+ + 7-) - 

[ire2 + 3n2A - ^£(27+ + 7-)r 

These expressions differ from the result (27) and the 
corresponding expression for F^o, For for the full- 
space reduced density matrix in the second or the first 
order of F, respectively, i.e., they both agree in the low- 
est order in F as should be expected. Thus, for small 
values of F we obtain very good agreement between 
the two approaches (curves are almost indistinguish- 
able from the solid lines in Fig. 2). This finding finally 
justifies the results obtained by the full approach as 
physically plausible although it does not yield any 
hints where the problem of the full approach might lie 
nor does it say anything reliable for larger F's. 



7. Discussion and conclusions 

Wo have presented our results for the dissipative 
DQD system and explicitly pointed out the paradoxes 
stemming from the Markovian GME description of 
this system, such as charge conservation breaking be- 
tween different junctions or unphysical negative values 
of the current noise. The weak coupling prescription is 
known to possess severe conceptual problems [2,17] in- 
cluding the ambiguity of the choice of the kernel (e.g., 
direct weak coupling prescription a la Bloch-Rodfield 
vs. RWA/secular approximation on top of that), 
breaking of positivity of the reduced density matrix 
within the Bloch-Redfield formalism, and breaking of 
the equations of motion (Ehrenfest theorem) within 
the RWA [18]. Thus, in some sense one shouldn't be 
surprised to find similar inconsistencies in the DQD 
study. Yet, there are also important differences be- 
tween the above mentioned effects and the present 
findings such as relative importance of the various 
discrepancies even for very small coupling constants, 
i.e., significantly increased sensitivity of these effects 



with the presence of the other bath. This indicates 

that the phenomena encountered here may be going 
beyond the "standard" weak coupling paradoxes and 
are specific to multiple-bath setups. 

This conjecture seems to be supported by further 
facts: the charge conservation is fulfilled without the 
presence of the dissipative bath which is a necessary 
condition for the occurrence of the breaking. RWA 
applied to our system expectedly breaks the equations 
of motion (physically incorrect symmetry between 
the emission and absorption processes), however, the 
full Bloch-Redfield kernel not only disobeys the posi- 
tivity (negative noise) but also breaks the equations 
of motion (charge conservation breaking) which does 
not seem to follow the standard weak coupling behav- 
ior. The multiple-bath dissipative systems have been 
known as challenging for quite long but at the same 
time they are characterized by physically interesting 
and sometimes counterintuitive behavior [19,20] which 
can be also responsible for the present findings. More- 
over, as a part of the diploma thesis of the first author 
[15] other systems were studied, e.g., the energy trans- 
port in two linearly coupled harmonic oscillators being 
an exactly solvable analog to the dissipative DQD but 
no similar phenomenon was observed there. This leads 
to speculations that these findings are determined not 
only by the simultaneous presence of more baths but 
also by the non-Gaussian character of the associated 
noise missing in the linear systems. 

To sum up, although the reported issue resembles the 
notorious difficulties and paradoxes of the weak cou- 
pling theory, neither its exact origin nor possible cures 
have been uniquely identified so far. Yet, we present our 
findings to the community in order to draw attention 
to this open nontrivial problem which may be, apart 
from its relevance for the particular physical system 
under study, of interest in the broader context of open 
dissipative quantum systems. In particular, our find- 
ings raise the questions of the development of charge- 
conserving approximation schemes within the gener- 
alized master equation approaches analogous to their 
non-equilibrium Green's function counterparts and of 
general understanding of dynamics of quantum sys- 
tems coupled to multiple baths. 
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